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Abstract: Weakly coupled, almost massless, spin particles have been pre- 
dicted by many extensions of the standard model of particle physics. Re- 
cently, the PVLAS group observed a rotation of polarization of electromag- 
netic waves in vacuum in the presence of transverse magnetic field. This 
phenomenon is best explained by the existence of a weakly coupled light 
pseudoscalar particle. However, the coupling required by this experiment is 
much larger than the conventional astrophysical limits. Here we consider 
a hypothetical self-interacting pseudoscalar particle which couples weakly 
with visible matter. Assuming that these pseudoscalars pervade the galaxy, 
we show that the solar limits on the pseudoscalar-photon coupling can be 
evaded. 
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1 Introduction 



In a recent paper Jain and Mandal [Tj proposed a mechanism to evade 
the astrophysical bounds on light pseudoscalars. The authors argued that 
there are no experimental or astrophysical bounds on the pseudoscalar self- 
couplings, which can, therefore, be treated as a free parameter. For suf- 
ficiently large, but perturbative, self-coupling the mean free path of pseu- 
doscalars can be very small inside the sun, provided we have sufficiently 
high density of these particles. The authors argued that the pseudoscalars 
produced inside the sun will accumulate due to self-couplings. The mecha- 
nism required the fragmentation of the pseudoscalars through the process 
4>{ki) + (j)(k 2 ) — > (j>(k 3 ) + 4>{ki) + 4>(k 5 ) + (j)(k 6 ). Some mechanism for energy 
loss as the particles propagate through the sun was also required in order 
that the pseudoscalar particles are at sufficiently low temperature. Using 
this mechanism the authors argued the results of the PVLAS experiment 
[2] can be consistent with the astrophysical bounds [3J and with the results 
of the CAST experiment [3]. The PVLAS collaboration finds a rotation of 
polarization of light in vacuum in the presence of a transverse magnetic field. 
If we interpret this rotation in terms of the coupling of a light pseudoscalar 
particle to photons we find that the allowed range of parameters are in con- 
flict with the astrophysical bounds [3J, although there is no conflict with 
any laboratory bounds [5]. The PVLAS result has motivated considerable 
theoretical work [6-15] as well as proposals for new experiments [T6l ITT] and 
observations [IB] . Furthermore new contraints have been imposed on ax- 
ion monopole-dipole coupling [19] assuming the mass range observed in the 
PVLAS experiment. 

Self-interacting dark matter has earlier been considered by many authors. 
It was proposed in Ref. [20] due to its interesting consequences for cosmol- 
ogy. Constraints on such a dark matter candidate and its implications for 
cosmology have been further studied [21J. The dark matter studied in the 
present paper may not give a significant contribution to the energy density 
of the universe due to its very small mass. Significant contribution may arise 
only if the pseudoscalar field undergoes coherent oscillations around the min- 
imum of the potential such as those expected for axion like fields [22]. This 
oscillatory field behaves just as non-relativistic dark matter. Some evidence 
for self-interacting dark matter is also found in considering the galactic dark 
matter distribution [23]. The standard cold dark matter scenario leads to 
cuspy dark matter galactic halos, which are not in agreement with observa- 
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tions [21]. The self interacting dark matter solves this problem and has been 
studied by many authors [25] . Alternatively in Ref . [26J it is claimed that the 
dark matter profile is in good agreement with cold dark matter predictions. 
It is interesting to determine whether the dark matter we consider can be 
the dominant component of the galactic dark matter. However we postpone 
this question for future research and for now simply assume that the pseu- 
doscalars we consider have densities negligible compared to the galactic dark 
matter density. 

In the present paper we propose a simple mechanism for evading the 
astrophysical bounds. Our basic assumptions are 

1. The pseudoscalars have a perturbative, but relatively large, self- coupling. 

2. The pseudoscalars form some component of the galactic dark matter. 

The second assumption implies that pseudoscalars at nearly zero temperature 
are present throughout the galaxy. This assumption is quite reasonable since 
pseudoscalars will be produced in the early universe and will be present 
today in the form of dark matter. Stars are formed in regions of dense 
molecular clouds where the average density of matter is much higher then 
the mean galactic density. It is very likely that the density of dark matter 
is also higher in these regions due to the large gravitational attraction. As 
the cloud of dust and gas collapses to form a star the pseudoscalar cloud 
will also collapse. This will lead to a high density of pseudoscalars inside 
the star. The mean free path of pseudoscalars produced inside the star may 
then be very small due to the their scattering on background pseudoscalars 
by the reaction (j){k\) + ^(fe) — > ^(^3) + <j>{k^). Hence this mechanism can 
trap pseudoscalars inside the star and and considerably limit the radiative 
transfer [27] that can occur through pseudoscalars. In the present paper we 
investigate this mechanism in order to determine whether this mechanism 
can indeed evade the astrophysical bounds. We will focus on the bounds 
imposed by considering the pseudoscalar production inside the core of the 
sun. Our mechanism also applies to bounds due to red giants and supernovae. 
However we postpone this discussion to future research. 
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2 Evading the Solar Bounds 



The effective Lagrangian for the pseudoscalars can be written as 

C = --F^F^ + -djd^ - -m 2 ct) 2 + -L-tFuvFt" - - x (f) A (1) 
4 2 M 2 4M^ M 4! v ; 

where we have assumed a self-coupling term besides their interaction with 
photons. Assuming that the results of PVLAS experiment can be explained 
in terms of pseudoscalar-photon mixing, the parameter and the pseu- 
doscalar mass m^, respectively, lie in the range 1 x 10 5 GeV < < 
6 x 10 5 GeV and 0.7 meV < < 2 meV [2]. We assume that the self- 
coupling a>\ = A 2 /47r < 0.1, so that it can be treated perturbatively. 

The density of dark matter in the galaxy may be parametrized as [28], 

P{r) = ~^~ 2 (2) 

We may assume that the pseudoscalars form a fraction £ of the galactic dark 
matter, where C Q = 4.6 x 10 8 M Q kpc~ 1 and a = 2.8 kpc. Using the mass 
of the pseudoscalar to be about 10~ 3 eV and the distance of earth from the 
galactic center r = 8 kpc, we estimate the number density n = 2.4£10 11 cm -3 
We assume £ to be less than 0.1 so that pseudoscalars contribute negligibly to 
the galactic dark matter density. This parametrization of the density profile 
is sufficiently reliable for our purpose. It correctly models the expected 1/r 2 
at large distances and becomes constant at small distances. The simulated 
dark matter distributions for several different models such as the standard 
cold dark matter and cold dark matter with self interaction are shown in 
Ref. [23] ■ The resulting density obtained for large r is in good agreement 
with that obtained with the parameterization given in eq. [2] for these models. 
Hence we find that a wide range of reasonable models of galactic dark matter 
produce a density in good agreement with what is obtained from eq. [2] at 
the position of the earth inside the milky way. 

The pseudoscalar density may be higher at the site of star formation due 
to the higher concentration of matter in this region. We set the pseudoscalar 
number density in region which lead to the formation of sun to be a factor rj 
times its galactic density. This region, which contained the primordial cloud 
of gas and dust whose collapse lead to the formation of the solar system, 
is expected to be of size larger than 100 - 1000 AU. As the star forms the 
pseudoscalars in this region may also collapse onto the star. The precise 
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pseudoscalar density profile will be determined later by requiring that the 
solar system, including the pseudoscalar gas, are in steady state. Here we 
next make an order of magnitude estimate of the pseudoscalar number density 
inside the sun, assuming that almost all the pseudoscalars within a radius 
Rs ~ 1000 AU collapse onto the sun. The resulting number density of 
pseudoscalars inside the sun is found to be 



w 2.4^10 



20 



1000AU 



cm' 3 (3) 



As the star collapses the temperature in the core increases, eventually 
leading to nuclear reactions. The photons in the core then start producing 
pseudoscalars through the Primakoff process pQ. The temperature of the 
pseudoscalars inside the sun also increases, partially due to the collapse of 
pseudoscalars and partially due to production of high energy pseudoscalars by 
the Primakoff process. We can compute the kinetic energy of pseudoscalars 
as they reach the solar radius from nearly infinite distance. Assuming that 
their initial kinetic energy is almost zero we find that the final kinetic energy 
is 3.4 x 10 -21 ergs which is equivalent to a temperature of 2.5 x 10 -5 K. 
Hence this predicts a relatively cool pseudoscalar gas. This result may be 
significantly modified due to conversion of photons into pseudoscalars inside 
the sun. As already mentioned the final density profile is obtained later 
assuming steady state of the solar system and does not rely on these rough 
estimates. 

The typical energy of a pseudoscalar produced inside the core of the sun, 
assuming present day conditions, is of order 1 KeV. This is much larger than 
the kinetic energy acquired by pseudoscalars due to conversion of potential 
energy. We expect that the pseudoscalars will eventually reach a steady state, 
where the outward pressure balances the gravitational attraction. Further- 
more there should be no net exchange of energy between the photon and the 
pseudoscalar gas. The rate at which photons loose energy by conversion to 
pseudoscalars is given by 

where p is the density of the medium and the angular brackets denote thermal 
averaging. Equating this to the rate at which photons are gaining energy due 
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to the conversion of pseudoscalars into photons we find 




do 



dE^ 



) 



P 




(5) 



It is clear that the maximum temperature that the pseudoscalars can acquire 
is equal to the photon temperature. The pseudoscalars would acquire this 
temperature if in thermal equilibrium with the photon gas. At the tempera- 
tures corresponding to the solar interior the pseudoscalars would be relativis- 
ts and their number density would be equal to the photon number density. 
This argument gives us a number density in the core of the sun to be of the 
order of 10 23 cm~ 3 . This is really a lower bound since, as discussed later, 
the temperature of pseudoscalars may be lower. In this case their number 
density has to be higher in order to maintain steady state with the photon 
gas. 



Figure 1: Feynman diagram for the fragmentation process 00 — > 0000. 

The pseudoscalars may not be in thermal equilibrium inside the sun. This 
is due to the fragmentation process 



which occurs through the Feynman diagram shown in Fig. 1. If we as- 
sume coupling A of order unity, the cross section for this process at center 
of mass energy E, a ~ a\/E 2 , is much larger than the cross section for 
the inverse Primakoff process for the conversion of pseudoscalars to photons. 
Hence pseudoscalars would fragment rather than convert to photons. The 
fragmentation process will stop only when the pseudoscalar energy becomes 




0(/ci) + 0(/c 2 ) -> 0(A; 3 ) + (f>(h) + Hh) + <f>(h) 
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comparable to their mass. Hence we expect the pseudoscalars inside the sun 
will have energies comparable to their mass and will not be ultra-relativistic. 

The cross-section o for the process <j>(ki) + 4>{ki) — ► 0(^3) + 0(^4) at 
leading order in perturbation theory is given by p] 



8EI 



1.2 x 10 n a A 



10" 6 GeV 

Errrt. 



GeV 



-2 



(6) 



Assuming a center of mass energy of 1 KeV, we find the mean free path 
/ m 2 x 10~ 7 /o;a cm. With ot\ = 0.1, this gives a mean free path I « 2 x 10 -6 
cm, which is much smaller than the mean free path of photons inside the core. 
Hence we expect that pseudoscalars will contribute negligibly to radiative 
transport in comparison to photons and this mechanism evades the limits 
imposed by considering cooling rate of the sun. 

As seen above, the pseudoscalars are essentially trapped inside the sun 
and unable to escape freely. Due to the fragmentation process they may 
not be in equilibrium with the rest of the system inside the sun. However 
in steady state there would be no exchange of energy between pseudoscalars 
and the remaining particles. Due to the fragmentation process a pseudoscalar 
of energy E will produce roughly E/m^ pseudoscalars. The typical energy 
of a pseudoscalar produced in the outer regions of the sun by conversion 
from photons is of order 1 eV. The temperature of the pseudoscalar gas is 
of order of their mass, i.e. of order 0.001 eV. Hence a photon produced by 
pseudoscalars will have energy of order 0.001 eV. Since the cross sections for 
Primakoff and the inverse Primakoff process are roughly equal, we expect 
that, in steady state, the number density of pseudoscalars in the outer re- 
gions of the sun is of order 1000 times the photon density. The pseudoscalars 
will, therefore, form a relatively cold and dense gas. For a certain range of 
self-coupling it cannot be described as an ideal gas since its range of interac- 
tion may be much larger than the typical separation between particles. The 
pseudoscalar gas will maintain steady state with the remaining particles by 
producing large number of low energy photons which will be quickly absorbed 
by the fermions inside the sun and hence will equilibrate. 

Outside the sun the photons undergo free streaming. The pseudoscalars, 
however, are unable to propagate freely due to the background distribution 
of pseudoscalars. We next get an estimate of the pseudoscalar density profile 
outside the sun by considering two limiting cases. We first assume that 
pseudoscalars are in exact thermal equilibrium with the remaining particles 



7 



inside the sun. In the other extreme we assume that the fragmentation 
process dominates. 



2.1 Thermal Equilibrium 

Here we assume that pseudoscalars are in thermal equilibrium inside the 
sun. This provides an important limiting case, although for the range of 
couplings that we are studying this may never be realized. In this limit the 
pseudoscalars are ultra-relativistic near the surface of the sun. The number 
density of pseudoscalars is obtained by eq. \5\ and is found to be same as that 
of photons throughout the radius of the sun. We can compute the mean free 
path of pseudoscalars throughout the solar radius using this density profile. 
It is found to be of order 10~ 7 /a\ cm near the core and increases to about 
10~ 4 /«a near the solar surface. Hence for a wide range of values of ot\ it is 
much smaller than the photon mean free path and the energy transport will 
occur primarily by pseudoscalar emission. Let R® denote the solar radius. 
For r > R Q , we assume, for simplicity, that we can describe pseudoscalars as 
an ideal Bose gas. At r = R Q we impose the boundary condition that the 
pseudoscalar temperature is equal to the temperature at the surface of the 
sun. The equation for the pressure gradient can be written as 



dP GM{r) , . n , 

= — ^rW' , ) + p «). (7) 

where is the energy density of the pseudoscalars, P the pressure and 

rr 

M{r) =Att dr'r' 2 P(f> (r') (8) 
Jo 

We point out that we are interested in the solution for r > R & , where we 
neglect all other particles except the pseudoscalars. In eq. [7] we have kept 
the dominant relativistic corrections. For r > R Q , we have P(r) = p^(r)/3 
and we may assume that M(r) — M & . The solution can be written as 



r) = p (i? )exp 



(9) 



This essentially means that p^r) remains roughly constant inside the solar 
system. The solution is not valid as we go to distances beyond the so- 
lar system since then we need to take into account gravitational effects of 
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other objects in the milky way. The pseudoscalar density profile obtained 
for r > R & again implies a very small mean free path in this region. Hence 
the pseudoscalars are unable to propagate freely throughout the entire solar 
system and the standard solar model remains unaffected by the presence of 
these particles. 

2.2 Fragmentation 

We next consider the other extreme where the fragmentation dominates. In 
this case we assume that the energy per particle of the pseudoscalars is of the 
order of their mass due to the fragmentation process. The pseudoscalars may, 
therefore, not be ultra-relativistic. In order to get an estimate of the density 
profile for r > R® we assume that they are non-relativisitic in the outer 
regions of the sun, such that their temperature is smaller than their mass. 
The temperature of the photon gas near the surface is about 10 4 K. Since the 
mass of the pseudoscalar is about 10~ 3 eV, we expect that the pseudoscalar 
temperature < 10 K. By using Eq. \5\ and the fact that the cross sections 
for the direct and the inverse Primakoff process are roughly equal, we find 
that the number density of pseudoscalars at the surface u^Rq) > 10 3 n 7 . The 
factor 10 3 arises due to the ratio T 7 /T^. We may similarly obtain an estimate 
of the density profile inside the entire radius of the sun. The precise value 
is not relevant for our purpose. The number density is certainly atleast as 
large as the photon number density throughout the sun and that is enough to 
suppress the energy loss through pseudoscalar emission to negligible values. 
The equation for hydrostatic equilibrium in this case is given by 



We next need to specify the nature of change that the medium undergoes. 
We shall be interested in the pseudoscalar gas for r > R Q . In this region 
we assume that it is governed entirely by its self-coupling and its interaction 
with other particles is negligible. Due to the large cross section of the process 
<p(f) —>■ cfxp it is clear that radiative transport of energy is negligible. This is 
because pseudoscalars are not able to propagate freely in this medium. The 
medium is therefore likely to undergo adiabatic changes and hence will reach 
a steady state such that Pp^ 1 = const, where 7 = cp/cy- For generality we 
assume a polytropic equation of state 



dP 

dr 



GM(r) P<t> (r) 



(10) 




(11) 
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which models a large number of processes for different values of the index 7'. 
For adiabatic change 7' = 7. We make the change of variables p^rp 7 '" 1 ' = 
z(r) and x = r/y/a where a = JCj' /[^(j' — 1)G]. In terms of x and z 
the equation for hydrostatic equilibrium can be brought into the form of 
Lane-Emden equation, 

hi + = ° < i2 > 

In order to bring the equation into dimensionless form it is convenient to use 
the scaled variables z and A, defined by, 

z(x) = z{x) P4> {R Q )^ 

x = \p^R Q )^- 2)/2 (13) 

where p^(R & ) is the pseudoscalar density at the surface of the sun. The 
equation in terms of z and A is same as eq. [12] with x and z replaced by A 
and z respectively. The boundary conditions to be imposed are 

z(K) = 1 
§<A.) = 6 (14) 

where A Q is the value of the dimensionless variable A at the surface of the 
sun r = R Q and 6 can be determined by relating dz/dX to dp^/dr at the 
surface of the sun. For the values given above, we find 5 ~ 1 x 10 -40 . The 
resulting density profile is plotted in Fig: [2j In obtaining this result we have 
set 7' = 7 = 5/3. In the present case the medium is essentially opaque to the 
propagation of pseudoscalars due to the large scattering cross section of the 
process <fxft —>■ 4><P- Hence we can safely assume negligible energy transport 
and use the adiabatic equation of state. 

We find that the density distribution is constant from the edge of the sun 
upto distance of order 100 Kpc. Ofcourse our solution is no longer valid once 
we are outside the solar system since we have neglected the presence of other 
objects in the milky way. Beyond a certain distance the pseudoscalar density 
has to match the galactic pseudoscalar density. However we do not address 
the issue of galactic density profile in this paper. We point out that even 
at very large distances from the sun, where the pseudoscalar density reaches 
its galactic values, the mean free path of pseudoscalars is still very small, of 
order of a cm, for a suitable choice of parameters a\ and £. Hence even here 
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Figure 2: Density Profile of Pseudoscalars With Distance. 

pseudoscalars may not be able to propagate freely for a wide range of these 
parameters. 

It is easy to see that the density has to remain constant for a large range 
of values of r above the surface of the sun. For this purpose we start with the 
basic equation for pressure gradient, eq. [TU1 Using the polytropic equation 
of state, eq. [HI we convert this into a differential equation for 

The total mass contained within a distance r, M(r), may be split up as 

rr 

M(r) = M Q + Anr' 2 dr' P4> (r') (16) 
Jr q 

The second term on the right hand side is negligible compared first term as 
long as r is comparable to or only a few orders of magnitude larger than 
Rq. We may, therefore, drop the second term for a wide range of values of 
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r. With this approximation the equation can be integrated analytically and 
we find 

Mrr'-*(^r' = -<7-l)^(£-i) (17) 

We find that the right hand side is negligible for r > R Q and hence for a 
wide range of values of r the solution is p^(r) ~ p^(R Q ). 

The number densities obtained for all values of r are again sufficiently 
large so that the mean free path of pseudoscalars is much smaller than that 
of photons. Hence the energy loss is expected to arise primarily through 
photon emission and the standard solar model will remain unaffected. 

It is interesting that the form of the solution for a wide range of values 
of r is same as that obtained in section 2.1 assuming thermal equilibrium 
throughout the sun. The number density is essentially constant within the 
entire solar system for r > R®. The precise number is obtained by the 
boundary condition imposed at the surface of the sun. Furthermore our 
mechanism for evading astrophysical bounds on pseudoscalar coupling and 
mass values is found to be applicable independent of the assumptions used 
to obtain the density profile inside the sun. 

We next check if our solution is consistent with the constraints imposed 
by the Pioneer anomaly (30]. We find that the solution obtained is not 
in conflict with the maximum allowed dark matter density inside the solar 
system [3T] [32] . However the density is about an order of magnitude small 
compared to what is required to reproduce the anomalous acceleration of 
pioneer [30J crafts 10 & 11, towards Sun. We may achieve higher densities if 
the pseudoscalars have temperature of the order of 1/2 K at the surface of the 
sun, instead of 10 K, which we assumed. This can arise if the pseudoscalars 
also loose energy as the propagate inside the sun [1] . By using Eq. [5] we find 
that this gives a higher pseudoscalar density at the surface of the sun by a 
factor of about 20. This leads to an increase in the density in such a way 
that it exactly matches the pioneer anomaly data. 

3 Medium modifications 

We next address an important consistency check of our proposal. It is impor- 
tant to determine how the presence of background pseudoscalars affect the 
results of the PVLAS experiment. The density of pseudoscalars is relatively 
small and due to their weak coupling with photons, they will not significantly 
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affect the propagation of photons. The dominant effect is modification of the 
pseudoscalar self energy. The medium modifications can be taken into ac- 
count by using the modified mass = m\ + Sm 2 h , where the thermal mass 
[33], 

5m* = / k 2 dk— -=± . (18) 

th 4tt 2 J E k e^ - 1 K ' 



where Ek = yk 2 + m 2 and (3 = 1/T. In the high temperature limit, i.e., 
T » m , the result is straightforward and it turns out to be 

8< = ^ (19) 
In the low temperature limit (3 — > oo, we find 



At the position of earth, we expect that the temperature would be at most 
as high as the mass of the pseudoscalar particle, due to the fragmentation 
process. Hence the change in the pseudoscalar mass due to thermal effects 
would be small as long as A < 1. The modification in mass could be large if 
for some reason the fragmentation process is not very effective. In this case 
we should interpret the mass extracted by the PVLAS experiment as being 
dominantly the thermal mass. In either case our mechanism for evading 
astrophysical bounds remains applicable. 



4 Summary and Conclusions 

We have shown that the solar astrophysical limits can be evaded provided 
we assume a sufficiently large value of the self-coupling of pseudoscalars and 
that pseudoscalars form some component of the galactic dark matter. The 
self-coupling required is within the perturbative regime and the galactic den- 
sity of pseudoscalars is assumed to be a negligible fraction of the dark matter 
density. Inside the sun the maximum possible value of pseudoscalar temper- 
ature at any radius r is equal to the solar temperature at that point. The 
corresponding mean free path of pseudoscalars is found to be sufficiently 
small so that they contribute negligibly to radiative transfer. We have shown 
that outside the sun the density profile of pseudoscalars is roughly uniform. 



13 



The density outside the star is found to be such that the mean free path 
of pseudoscalars is smaller than a cm for a wide range of values of the self 
coupling a\. Hence the pseudoscalars cannot escape freely from the sun 
and the standard solar model remains unaffected. The solar limits on the 
pseudoscalar-photon coupling are, therefore, evaded. 
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